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Summary 

Simulations with Aarseth's (1994) NB0DY5 code are presented of an initially dense binary-rich 
cluster. It is assumed that the star formation efficiency is 50 per cent with instantaneous mass 
loss. 

The model central density and velocity dispersion agree with the observational constraints if 
expansion is only about 6 x 10^ yr old. Additionally, the observed binary proportion constrains 
the primordial proportion to have been significantly less than in Taurus-Auriga. 

The claim that a variation of the birth binary proportion with gas cloud parameters 
has been detected, however, can only be verified if the cluster can be shown to be expanding 
rapidly. 



^Poster presented at the NATO ASI on The Physics of Star Formation and Early Stellar Evolution (C.Lada & N.Kylafis), 
Crete, 24.05.-5.06.98. Further details and additional models will appear in Kroupa, Petr & McCaughrean (1998). 
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1. Introduction 

The Trapezium Cluster in Orion is an important astrophysical laboratory for the study of the formation of 
stars and planetary systems, and of the early dynamical evolution of a young cluster rich in binary systems. 

• It's age is estimated to be younger than 5 Myr, and most probably only a few 10^ yr. 

• The central stellar number density is unusually large for embedded clusters, with mean 
inter-stellar separations of about 6000 AU. 

The distribution of orbital elements of a population of binary systems carries a memory of the 
dynamical history of this population. Short-period primordial binary systems are an energy reservoir, 
whereas long-period binary systems are an energy sink. The cluster and its constituent binary systems 
interact, which in part drives cluster evolution. 



• Surveys of Galactic field main-sequence steirs show that a significant fraction (/tot ~ 50%) are in 
binary systems (e.g. Duquennoy & Mayor 1991). 

• Observations of sparse groups of pre-main sequence stars in Taurus- Auriga imply a significant 
overabundance (/tot ~ 1) of binary systems (e.g. Kohler & Leinert 1998). 

• In the Trapezium Cluster the binary proportion is comparable to the Galactic field value (Petr et 
al. 1998a,b, Prosser et al. 1994). 



QUESTION: Can a Taurus-Auriga binary proportion have dynamically evolved to the reduced value, 
or is the Trapezium population primordial ? 
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2. Models 
2.1. The initial model Trapezium Cluster 

• N = 1600 stars (point masses) (McCaughrean & Stauffer 1994). 

• Half-mass radius i?o.5 = 0.1 pc (McCaughrean & Stauffer 1994). 

• IMF from Kroupa et al. (1993) for m < 1 Mq, and Scalo (1986) for m > IMg. 

• Lower and upper stellar mass limits of m\ = 0.08 M© and mu = 30 Mq, respectively. 

• Plummer density distribution. 

• Initial position and velocity vectors are independent of stellar mass. 

• Velocity distribution of the binary centre-of-masses is isotropic. 

The resulting cluster mass is Md = 700 M©. If the Trapezium Cluster were in virial equilibrium then the 
model relaxation and crossing times would be treiax = 0.62 Myr and fcross = 0.1 Myr, respectively. 

The virial ratio Q = ii^kin/l^'potl = 0.5 in virial equilibrium. Here Q = 1. 

The initial velocities are chosen to correspond to a system with a combined mass 
Mstars + Mgas = 2 X Mgtars- That is, it is assumed that the massive stars in the cluster have 
driven out a gas mass, Mgas ~ Afgtars = Mc\ = 700 Af© (cf Churchwell 1997), immediately after they "turn 
on" and before any significant stellar-dynamical processes occur. This corresponds to a star-formation 
efficiency of 50 per cent. 
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2.2. Primordial binary systems 

The total (summed over all periods) binary proportion is 

^tot — ,r ,T I l-^j 

JVbin + JVsing 

where iVbin and Nsiag are the numbers of bound binary and single star systems, respectively. 

• /tot — 1 assumes that the binary-star properties do not vary with star-forming conditions, 

apart from the effects of crowding, and that they are identical to what is observed in Taurus-Auriga. 

• /tot = 0.6 assumes that the binary-star properties vary with star-forming conditions in the 

sense suggested by Durisen & Sterzik (1994). 

The initial model period distributions are compared with the observational data in Fig. 

Additionally, the following assumptions are made: 

• The initial mass-ratio distribution is obtained by random pairing of the stars. 

• The initial eccentricity distribution is thermally relaxed. 



Fig. |2| shows the resulting initial eccentricity-period diagram, after pre-main sequence eigenevolution. 
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3. Results 

• The central density and tangential velocity dispersion decrease rapidly due to the expansion. 
They are consistent with the observed values at time t « 0.06 Myr. 

• Rapid expansion halts the disruption of binary systems at an early stage. 
This is shown in Figs. p| and ^ 
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Implications: 

• The expanding model could be a reasonable description of reality if about 50 per cent of the mass of 
the cluster was expelled about 60 thousand years ago, at which time the present Trapezium Cluster 
would have gone into a rapid expansion phase. This solution implies an extreme youth of the 
gas-free Trapezium Cluster. 

This would be consistent with claims that the circum-stellar material seen around the young stars 
in the Trapezium Cluster should be removed within 10^ — 10^ yr through photo-ionization by the 
most massive central star, 9^ C Ori (e.g. Bally et al. 1998). 

• The model apparent binary proportion is consistent with the observational constraint for 
t ^ 0.05 Myr, provided /tot ~ 0.6 initially! 

This is lower than the binary proportion in Taurus— Auriga. If model B2 does represent reality, 
then this may be due to dynamical evolution in the proto-cluster prior to gas expulsion, or due to a 
dependency on cloud temperature as suggested by Durisen & Sterzik (1994). 

Main shortcoming: 

In the present study the Trapezium Cluster is treated as an isolated entity. 

In reality, it appears to be the core of the much more massive and extended ONC, which is partially 
embedded in the parent elongated molecular gas cloud (Hillenbrand & Hartmann 1998). 
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Table 1: Initial conditions for the Trapezium Cluster models. Three simulations are performed 
for each model. Column 5 contains the mEiximum binary-star period. The minimum period is 

-Pmin = 1 day in all cases. The form of the period distribution is defined in column 6, where K2 (Kroupa 
1995a) refers to a birth distribution that is consistent with young systems in Taurus- Auriga, and DM91 
(Duquennoy & Mayor 1991) refers to the Gaussian log-period distribution of the Galactic field. A 1 Mq 
system has, with logioPmax = 8.43 (11) a semi-major axis of 8200 AU (4.2 x 10^ AU). 
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Fig. 1. — Distribution of orbits, fp, for main sequence multiple systems (solid dots, Duquennoy & Mayor 

1991) and prc-main sequence systems in Taurus- Auriga (open squares; logio-P > 4: Kohler & Leinert 1998, 
logio-P = 3.5: Richichi et al. 1994, logioP < 2: Mathieu 1994). The dotted histogram is the initial period 
distribution from Kroupa (1995a). Crowding in the model Trapezium Cluster changes this distribution to 
the long-dashed one. A Gaussian log-period birth distribution that fits the solid dots, changes in the model 
Trapezium Cluster due to crowding to the distribution shown as the short-dashed line. 
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Fig. 2. — The initial eccentricity— period diagram for /tot = 1- The distribution for /tot = 0.6 is the 
same in eccentricity, but has fewer orbits for fogioi-" > 4. The thick dashed Une represents the observed 
envelope for main-sequence binary stars with a G-dwarf primary (Duquennoy & Mayor 1991). 
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Fig. 3. — Upper panel: The time evolution of the number of stellar systems within a spherical 
distance of i? = 0.053 pc of the position of the density maximum of the model chist(T. The thick curves 
assume no binary systems are resolved, and the thin curves count all stars. Lower panel: The time 
evolution of the velocity dispersion of ccntre-of-masses within a projected distance of r = 0.41 pc of 
the position of the density maximum of the model cluster. In both panels, the solid curves are for initially 
/tot = 1 (model Bl), and the dot-dashed curves are for initially /tot = 0.6 (model B2), and the observational 
constraints with the Poisson error range are indicated by the dotted lines. 
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Fig. 4. The time evolution of the apparent binary proportion in the model Trapezium Clusters. 
The solid line is for initially /tot = 1 (model Bl) and the dot-dashed line is for initially /tot = 0.6 (model B2). 
Upper panel: observational constraints from Prosser et al. (1994, r = 0.249 pc, di — 26 AU, d2 — 440 AU) 
are shown as dotted lines. Lower panel: observational constraints from Petr et al. (1998a, r = 0.04 pc, 
di = 63 AU, d2 = 225 AU) are shown as the horizontal lines. Here, the central dotted line is /app for 
all systems in their sample, and the central dashed line is /app for the low-mass systems only. Poisson 
uncertainties are indicated by the upper and lower horizontal lines. 



